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Abstract 



We compare predictions for the spin-independent contribution to the neutralino-proton 
scattering cross section a xp and for the anomalous magnetic dipole moment of the muon, 
a n = (<?u — 2)/2, in models with gravity-mediated supersymmetry breaking. We nearly 
always find a positive correlation between these two measurables, i.e. scenarios with 
larger a M also tend to have larger a XP , but the detailed prediction differs greatly between 
,^Zh • models. In particular, we find that for the popular mSUGRA scenario with universal 

soft breaking masses at the scale of Grand Unification, measurements of a p currently 
seem more promising. On the other hand, if scalar soft breaking masses at the GUT 
scale receive sizable contributions from 50(10) D-terms, one often finds scenarios with 
large a xp but a M below the currently foreseen sensitivity. A string-inspired model with 
non-universal scalar spectrum at the GUT scale falls between these two cases. 



1 Introduction 



The Minimal Supersymmetric Standard Model (MSSM) |IJ is one of the best motivated ex- 
tensions of the Standard Model. It offers a natural solution of the hierarchy problem as 
well as amazing gauge coupling unification [H]. Since naturalness requires that at least some 
superparticles have masses at or below the TeV scale, supersymmetric theories generally pre- 
dict a rich phenomenology at future colliders [|J. As extra "bonus", the simplest version of 
the MSSM, where R— parity is conserved, also contains a new stable particle (the lightest su- 
persymmetric particle, LSP); in most cases this is the lightest neutralino, which often makes 
a good Dark Matter candidate ||. 

Unfortunately, direct searches for superparticles at high energy colliders so far yielded null 
results. Although it is quite possible that direct evidence for supersymmetry will be found at 
the upcoming "run 2" of the Tevatron collider, a decisive test of weak scale supersymmetry 
may only be possible at the LHC, which will not commence operations for another five or six 
years. In the meantime it is important to explore slightly less direct methods that might either 
give evidence for supersymmetry, or else constrain the (usually quite large) parameter space 
of supersymmetric models. Here we study two kinds of experiments which are expected to 
yield results with greatly improved sensitivity in the next couple of years: direct searches for 
ambient Dark Matter^], and the measurement of the anomalous dipole moment of the muon. 

Current direct Dark Matter search experiments are sensitive to scenarios with X\P scat- 
tering cross section a xv of order 10 -6 pb or more. In fact, the DAMA experiment even claims 
evidence for WIMP Dark Matter, based on the annual modulation of their counting rate ||; 
however, the CDMS experiment, which has comparable sensitivity, sees no signal 0. Within 
the next two or three years we can expect results from the upgraded CDMS experiment (at the 
deep underground Soudan site) || as well as from the CRESST experiment in the Gran Sasso 
laboratory || . The sensitivity of these experiments should be about two orders of magnitude 
better than that of current experiments, so that scenarios with a xp exceeding 10~ 8 pb might 
be testable. 

SUSY models can also be tested using precision measurements at low energy experiments. 
Through their loop effects, sparticles contribute to low energy physics and these effects may 
become significant if the masses of sparticles are not too large. The supersymmetric contribu- 
tion to the muon magnetic dipole moment (MDM) a p = \{g — 2) M is one of the most robust 
probes. Unlike supersymmetric contributions to rare decays or CP violation in the B— system, 
which are currently being studied at the B— factories BaBar and BELLE, the supersymmetric 
contribution to is not very sensitive to the way the supersymmetric flavor and CP problems 
are solved. At present, the muon MDM is measured to be a e * p = (116 592 05 ±45) x 1CT 10 [pLOf] 
and is consistent with the Standard Model. The ongoing analysis of data from the Brookhaven 
experiment E821 is expected to reduce the uncertainty further and the ultimate goal of the 
experiment is Aa M ~ 4 x 1CT 10 [[IlJ.Q Supersymmetric contributions can easily exceed this 



sensitivity without violating any other constraints on the supersymmetric particle spectrum 



Many physicists would accept direct detection of weakly interacting massive particles (WIMPs) as (com- 
ponent of) the dark halo of our galaxy as fairly direct evidence for weak-scale supersymmetry. However, 
given that direct detection experiments can basically only measure the mass and scattering cross section of 
the WIMP, such a detection, while undoubtedly of the greatest importance, may not convince everybody that 
supersymmetry has been found. 

2 In order to fully exploit such a small experimental error it is also necessary to reduce the hadronic un- 
certainty of the SM prediction, which currently stands at ~ 7 x 1CP 10 |HJ Here the analysis of data 
that is being collected at the low-energy e + e _ colliders in Frascati, Novosibirsk and Beijing, as well as from 
semi-leptonic decays of the r lepton, will be crucial. 
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(from direct searches or other loop effects) ||13|| . 

These supersymmetric contributions to stem from smuon-neutralino and sneutrino- 
chargino loops. Their size thus depends on the parameters that appear in the chargino and 
neutralino mass matrix, as well as on the masses of the sleptons. Of course, the parameters 
of the neutralino mass matrix also determine the mass and composition of the lightest neu- 
tralino, which we assume to form the Dark Matter in the Universe. Moreover, at least for the 
theoretically most appealing case of a bino-like LSP, the LSP relic density is essentially de- 
termined by the mass of the LSP as well as the masses of the SU(2) singlet sleptons |TJ]]. The 
predictions for the relic density of bino-like LSPs and for the supersymmetric contribution to 
a M are thus related. 

There also exist relations between the SUSY contribution to a M and the predicted LSP- 
nucleon scattering cross section. As already mentioned, both quantities depend on the param- 
eters of the neutralino mass matrix. An even more direct connection comes from the common 
dependence on the ratio tan/3 of the vevs of the two neutral Higgs fields of the MSSM, since 
both and a xp increase with increasing tan/3. The operator in the effective Lagrangian that 
gives rise to the anomalous magnetic moment of the muon couples left- and right-handed 
muons, i.e. violates chirality. In the MSSM (as well as in the SM) all chirality violation in 
the (s)muon sector is proportional to the Yukawa coupling of the muon. This coupling in turn 
is proportional to 1/ cos/3, which scales like tan/? if tan 2 /? ^> 1. The leading contributions 
to the spin-independent (coherent) contribution to a xv also involve violation of chirality (15 



but this time in the (s) quark sector .0 These contributions either come from the exchange of 
CP-even Higgs bosons, or from squark exchange. In both cases chirality violation ultimately 
comes from the quark Yukawa coupling, which, in case of down-type quarks, has the same 
tan/? dependence as the Yukawa coupling of the muon. Note that searches for Higgs bosons at 
LEP already imply tan/3 > 2 in the MSSMfJj; this implies that the Yukawa couplings of up-type 
quarks, which scale oc 1/ sin f3, must be close to their SM-values, and essentially independent 



of tan/3. The upshot of this discussion is that both the SUSY contribution to a M [OL W\ and 
°xp EH' 0' can ^ e l^te large if tan/3 ^> 1. 

In this paper, we investigate these connections between the SUSY contribution Aa /t to 
the magnetic dipole moment of the muon and the spin-independent neutralino-proton cross 
section quantitatively. We study the minimal supergravity (mSUGRA) model as well as some 
other SUSY models where the assumption of strict scalar mass universality at the GUT scale 
is relaxed. In Sec. 2, we discuss Aa^ and a xp in the mSUGRA scenario. These observables 
are below the sensitivity of near-future experiments in the small tan /3 region, partly because 
the experimental Higgs mass bounds then tightly restrict the allowed SUSY parameter space. 
However, they can become significant, and could be detectable, if tan/3 is large. In Sec. 3 we 
relax the assumption of universal scalar masses. We find that the correlation between Aa^ 
and a XP quite sensitively depends on the boundary condition at the GUT scale. Sec. 4 is 
devoted to conclusions. 



3 There are also squark exchange contributions where chirality violation is provided by the LSP mass, but 
these contributions are suppressed by an extra power of (rn~ — m|) _1 compared to the leading squark exchange 
contribution, and can thus usually be ignored in the type of model we are considering |jl5| . 

4 Higgs searches also permit tan/3 ^ 0.6, but this is strongly disfavored on theoretical grounds, since the top 
Yukawa coupling would have a Landau pole quite close to the weak scale. 
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2 mSUGRA 



In the minimal supergravity model, it is usually assumed that all squared scalar masses receive 
a common soft SUSY breaking contribution itl 2 q = mfj = m 2 D = m\ = m 2 E = rn 2 Hu = w^ H = 
m 2 at the GUT scale Mx — 2 • 10 16 GeV, while all gauginos receive a common mass M and 
all trilinear soft terms unify to A. The renormalization group (RG) evolution of soft breaking 
squared Higgs masses then leads to consistent breaking of the electroweak symmetry, provided 
the higgsino mass parameter \i can be tuned independently |T9|] . In this paper, we chose the 
weak scale input parameters m^im^) = 4.2 GeV, m t (m t ) = 165 GeV, and tan/3. We minimize 
the tree level potential at renormalization scale Q = ^/mfmt, which essentially reproduces the 
correct value of \i obtained by minimizing the full 1-loop effective potential [pOfl . 

With these assumptions, the mSUGRA model allows four continuous free parameters 
(m, M, A and tan/?). Also, the sign of \x remains undetermined. In Fig. 1, we plot the 
SUSY contribution to the muon MDM, Aa M , vs. the spin-independent neutralino-proton 
cross section, a xp , for three different choices of tan/3 = 4, 10 and 30. Here, we take A = and 
\i > and allow m and M to vary in the intervals m < 500 GeV and M < 600 GeV. We 
include loop corrections to the masses of neutral Higgs bosons from third generation quarks 



and squarks, including leading two-loop corrections [21]. We use the expressions of ref. 



for the calculation of Aa M . The calculation of a xp is based on refs. [fq, |22H. We use the value 



m s (p\ss\p) = 130 MeV for the strange quark's contribution to the nucleon mass; this matrix 
element is uncertain to about a factor of 2, leading to a similar uncertainty in the prediction 
of a XP . Finally, the calculation of the scaled LSP relic density fi x h 2 uses results of refs.jn], |23| ; 
s— channel poles are treated as described in ref.pi|. The co-annihilation of Xi with sleptons 
is not included. This effect can increase the cosmologically allowed region of parameter space 
towards higher LSP masses, if m^o ~ mf however, in this (limited) region of parameter space 
both Aa M and well below present and near-future sensitivity. 




°xp (pb) 



Figure 1: Aa M vs. a xp in mSUGRA for tan/3 = 4, 10,30 (from bottom to top region). We take A = and 
/i > 0, and scan m G [0, 500] GeV and M € [0, 600] GeV, subject to experimental constraints. The heavily 
marked points satisfy the requirement Q, x h 2 < 0.3. 
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Accelerator bounds significantly limit the SUSY parameter space. For the parameter range 
scanned in this plot we find that the lightest Higgs boson h couples essentially like the single 
Higgs boson of the SM; we thus demand that its mass rrih > 111 GeV. This follows from recent 



LEP results p5j, allowing for a 2 GeV theoretical uncertainty in the calculation of pi 



We further require that the chargino mass m-+ > 100 GeV |[26|| . We exclude regions where 
the LSP is charged, i.e. where < m x o or < m%o. The heavily marked points satisfy 
the further requirement Q x h 2 < 0.3. This bound on the matter density in the Universe follows 
from the analysis of recent cosmological data p7| . 



The lowest dotted region in Fig. 1 corresponds to tan/3 = 4. In this case both Aa M and a xp 
are too small to be detected in the near future. Furthermore, the lower bound on the Higgs 
mass forces one into a region of parameter space where the dark matter density is quite large 
and therefore cosmologically disfavored. Through radiative corrections, which increase with 
increasing stop masses, rrih depends quite strongly on M, but only weakly on m. The reason 
is that contributions from Yukawa interactions to the RG equations reduce scalar masses at 
the weak scale as compared to their GUT scale values. As a result, in the expressions for the 
weak scale squared stop masses and m 2 , m 2 appears with coefficients significantly less 
than 1, while M 2 appears with coefficients well above 1 [|l^]. For tan/3 = 4, the Higgs mass 
bound rrih > 111 GeV requires M > 500 GeV, which corresponds to m^o > 217 GeV. Such a 
large SUSY mass scale gives too large a DM mass density 



As tan/3 increases, both a M and a xp increase. As mentioned above there are essentially 
two types of diagrams involving superparticles which contribute to a^, i.e., neutralino-smuon 
and chargino-sneutrino loop diagrams. Since a M requires chirality violation, for tan/3 3> 1 
the dominant contributions are proportional to the product of an electroweak gauge coupling 
and the Yukawa coupling of the muon, where the latter factor either comes directly from the 
higgsino component of the chargino or neutralino in the loop, or from //£ — fiR mixing. For 
tan/3 ^> 1 one thus has Aa M oc tan/3. If mass splittings between different sparticles are not 
too large, so that their mass scale can be described by the single parameter msusY, the total 
result can be estimated as |16| 

l Ac Vl = ^"2 (aSH + a9i) ^^tan/3 (1) 
61tx a Kb o / mg USY 

where g± and #2 are the U{1)y and SU{2) gauge couplings, respectively. In our convention, 
where the gaugino masses and tan/3 are positive, the sign of Aa M is equal to the sign of fi. 

On the other hand, three classes of diagrams contribute to LSP-quark scattering: the 
exchange of a Z or CP-even Higgs boson in the t channel and squark exchange in the s 
or u channel. LSP interactions with matter can naturally be separated into spin-dependent 
and spin-independent parts. Z and squark exchange contribute to the former, and squark 
and Higgs exchange to the latter. In most situations the dominant contribution to the spin 
independent amplitude is the exchange of the two neutral CP-even Higgs bosons, although in 



some cases the contribution of squark exchange is substantial [13]. For given inputs m,M,A 
and sign(/i), the mass of the heavy CP-even Higgs boson decreases as tan/3 increases. This 
effect originates from the contribution of the bottom Yukawa coupling to the RGE of m 2 ^ |28 



recall that this Yukawa coupling also scales like tan/? for tan/? ^> 1. Moreover, its coupling to 
d— type quarks is essentially proportional to tan/3 for tan 2 /3 ^> 1 and m\ > m\ max , where 
is the mass of the CP-odd Higgs boson and m\ max is the rriA — * oo limit of the mass of the light 
CP-even Higgs boson. Under these conditions, which are almost always satisfied in the models 
we are studying, h couples essentially like the single Higgs boson of the SM, independent of the 
value of tan/?; rrih also becomes almost independent of tan/3 for tan/3 > 10. The contribution 
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to a xp from H exchange, which grows quickly with increasing tan/?, thus starts to dominate 
over that from h exchange once tan/? > {m H /m h ) 2 . Moreover, the leading down-type squark 
exchange contribution to a XP , which requires chirality breaking in the (s)quark sector, also 
increases oc tan/3 for large tan/3. 

We can see these enhancements in Fig. 1. For tan/3 = 10, Aa M is, in the cosmologically 
favored region, well above 1 x 10~ 9 and can reach upto ~3x 10~ 9 , which could be detected 
in the near future. Also, a xp increases quite a lot, but still remains below 10~ 8 pb. These 
enhancements, compared to the tan/3 = 4 case, partly come from the lower SUSY mass 
scale that is permitted by the Higgs search limit. For tan/3 = 10, the Higgs mass bound 
rrih > HlGeV only implies M > 260 GeV, which corresponds to m x o > 110 GeV; this allows 
scenarios with cosmologically favored DM mass density in the low m and M region. For 
tan/3 = 30, in the cosmologically favored region Aa^ takes values of 4 ~ 9 x 10~ 9 , while a XP is 
mostly above 1 x 10~ 8 pb and can reach upto 5 x 10 -8 pb. The smallest allowed SUSY mass 
scale for tan/3 = 30 is only slightly lower than that for tan/3 = 10, since, as stated above, 
becomes quite insensitive to tan/3 for tan/3 Z 10. 

Note that the highest Aa^ values shown in Fig. 1 already reach the current level of sensi- 
tivity, whereas the maximal value of a xp in this figure is still more than an order of magnitude 
below current sensitivity. We conclude that, in view of near future detect ability, Aa M is more 
sensitive than a xp in the mSUGRA model; this agrees with results of ref. |29|] (for m < 500 
GeV). One of the important factors for determining the size of a xp is the higgsino component 
of the LSP. MSUGRA predicts a bino-like LSP Xi f° r moderate values of m and M (below 
~ 500 GeV) |30|, [T4|]. This is a rather model independent result pi[ . Large positive cor- 
rections to squark masses from gaugino loops, together with the large top Yukawa coupling, 
drive the squared soft breaking Higgs mass m 2 Hu negative at the weak scale. On the other 
hand, correct symmetry breaking requires rn 2 H + /x 2 > — M|/2. One has to make large in 
order to obtain the correct electroweak symmetry breaking scale, if scalar masses and gaugino 
masses are of the same order. Since the Higgs-LSP-LSP couplings require higgsino-gaugino 
mixing, they scale like 1/ ' \i for fi 2 3> M§ . The Higgs exchange contribution to a xp then scales 
like 1/yU 2 . Similarly, the prediction for Aa^ depends on \i through gaugino-higgsino mixing 
as well as [l^ — jin mixing. The results of Fig. 1 thus depend quite sensitively on the value 
of fi at the weak scale, which in turn depends on the GUT scale boundary conditions. The 
correlation between Aa^ and a xp could also be different in generalizations of mSUGRA. That 
is the subject of the next section. 

Before going to the next section, we wish to discuss the dependence of our results on the 
parameter A as well as the sign of /i. The A— dependence is mostly indirect. First of all, A has 
some impact on rrih, and hence on the allowed region of the (m, M) plane. This can be seen 
in Fig. 2 a), where the Higgs mass dependence on M is plotted for three different choices of 
A = +2M, 0, — 2M. For a given M, a larger A gives larger m^. Therefore, if we chose a large 
positive A value, the minimal allowed value of M will be lowered, so that the maximal allowed 
value of Aa^ will be raised, even though this observable is almost independent of A if the 
slepton and -ino mass parameters are held fixed. A also has some impact on the [i parameter. 
As we can see in Fig. 2 b), larger A values give higher fi values and therefore more bino-like 
LSP for a given M. This reduction of the higgsino component of the LSP would reduce the 
Higgs exchange contribution to a xp . Some numerical examples of the effects of varying A are 
given in the next section. 

Changing the sign of [i from positive to negative has both indirect and direct effects. Since 
now A and /i appear with opposite sign in the t L — Ir mixing terms, the predicted value of 
rrih is reduced slightly. One thus has to increase the minimal allowed value of M in order 
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a) b) 

Figure 2: The dependence of a) m h and b) /i on M for three different choices of A = +2M, 0, — 2M (from 
top to bottom). Here we fix m=100 GeV for definiteness. The dotted and solid lines are for tan/3 = 5, 10 
respectively. 



to satisfy the bound nih > 111 GeV. In fact, we didn't find any allowed solutions within the 
scanning range for tan/? = 4, Aq = and /i < 0, while the minimal allowed value of M for 
tan/3 = 10 is increased by ~ 40 GeV, leading to a reduction of the maximal possible |Aa M | 
to about 2.4 x 10~ 9 . Note that for fixed values of the other parameters, the magnitude of 
Aa M is almost independent of the sign of //, but the sign of Aa^ changes when the sign of \i 
is flipped. In contrast, the sizes of the XiXi(h,H) couplings do depend quite significantly on 
the sign of /x, unless tan/3 3> 1 fl4j| . In particular, for /j < strong cancellations occur [|l^ 



both within different contributions to the same coupling, and between the h and H exchange 
contributions to a XP . As a result, for tan/? = 10, cross sections as low as 4 x 10~ n pb are 
possible in the cosmologically favored region, and the maximal allowed value of a xp is reduced 
by more than one order of magnitude, to 5 x 10~ 10 pb. The importance of the sign of /i 
diminishes with increasing tan/?; hence results for tan/? = 30 and \x < are quite similar to 
the results for tan/3 = 30 shown in Fig. 1. Finally, we mention that in models with universal 
scalar masses, solutions with fi < tend to have too large a branching ratio for radiative 
b — ► 57 decays |3^|.Q However, this prediction is sensitive to details of the flavor structure of 
the soft breaking terms, unlike the quantities we have depicted in Fig. 1. We therefore do not 
attempt to analyze the constraint from b — ► 57 decays quantitatively. 



5 Recently the leading higher order corrections to this branching ratio have been computed |Q. While 
a detailed parameter scan has not yet been performed, the published results indicate that scenarios with 
degenerate scalar masses and /i < remain problematic at large tan/3. On the other hand, for [i > agreement 
with the data can be obtained even for tan/3 ~ 45. 
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3 More general models 



In this section we relax our assumptions, allowing for non-universal soft scalar masses at 
the GUT scale, while keeping the unification of the gaugino masses. As specific models, we 
consider a string-inspired supergravity model and an 50(10) Grand Unified model. 

Soft SUSY breaking terms have been studied within the framework of string-inspired 



supergravity models pi. In a wide class of models one obtains the following relations among 



soft SUSY breaking terms: 

A ijk = —M, m 2 + m] + m 2 k = M 2 , (2) 

for non- vanishing Yukawa couplings Y^-/. of chiral superfields $ J and These relations 
hold if the theory is target-space duality-invariant and the Yukawa couplings in the supergrav- 
ity basis are field-independent. In addition we assume that only the F— terms of the dilaton 
and the moduli superfields contribute to SUSY breaking, and that the vacuum energy van- 
ishes. The generic form of supergravity theories leading to these relations has been obtained 
in Ref. These relations allow non-universality of soft scalar masses, but their overall 

magnitudes are bounded by the gaugino mass through the sum. This point is significant. 

Here we apply the relations (U) to the up and down sectors of (s) quarks and the (s)lepton 
sector separately. To be explicit, we have the relations between the soft breaking contributions 
to the squared scalar masses, 

m 2 Q + ml + m 2 Hu = M 2 
m 2 Q + m 2 D + m 2 Hd = M 2 

m 2 L + m 2 E + m 2 Hd =M 2 , (3) 

and a universal A-term, A = —M. 

In fact, these relations include quite a large parameter space, where sfermion masses are 
not equal to each other. However, as the first trial we take universal sfermion masses; this 
remains the simplest solution of the supersymmetric flavor problem. We shall later comment 
on a case with non-universal sfermion masses. 

The 50 (10) theory incorporates a complete generation of MSSM matter superfields into 
the 16-dimensional spinor representation, \1/ 16 . In addition to these matter superfield, the 
minimal 50(10) model includes a 10-dimensional Higgs superfield $io which contains the two 
Higgs superfields of the MSSM (as well as their SU(3) triplet, SU{2) singlet partners). When 
5O(10) breaks to the MSSM gauge group SU(3) C X SU{2) L x U{1) Y , additional D-term 
contributions (parameterized by M|, which can be either positive or negative) to the soft 
SUSY breaking masses arise p6| : 

m Q = m \ = m u = ™i6 + M D 
m 2 D = m 2 L = m\ % - 3M 2 D , m\ d = m 2 10 =F 2M 2 D , (4) 

where mi6 and mio are scalar soft breaking masses for fields in the 16 and 10 dimensional 
representations of 5O(10), respectively. 

The modifications (B|) and (Q) of the mSUGRA boundary conditions change our predictions 
in two different ways: through modifications of the slepton spectrum (which change Aa M and 
Q x h 2 ), and through the changed value of at the weak scale, which mostly affects a xp . While 

6 These relations are RG-invariant at one-loop for a single gauge group. 
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the changes of the slepton spectrum can be read off directly from eqs.(^) and the changes 
of are more subtle. The main effects can be understood as follows. In the mSUGRA 
scenario, the contributions of m and M to the weak scale values of the soft breaking Higgs 
boson masses can be parameterized as 

m 2 Hd ~ m 2 + 0.5M 2 ; 
m 2 Hu ~e H m 2 -3.3M\ (5) 

where we have assumed sin f3 ~ 1 but ignored contributions from the bottom Yukawa coupling. 
The coefficient e# is small, because the GUT scale value of m 2 H is canceled by the scalar masses 
appearing in the RG running.[] The effect of non-universality on m 2 Hu can be parameterized 
by introducing m 2 Q + mfj + m 2 Hu = 3m 2 and 5m 2 H = m 2 Hu — m 2 : 

m 2 Hu = 5m 2 H + e H m 2 s - 3.3M 2 . (6) 

This equation follows because the radiative correction to ra 2 Hu is proportional to m 2 Q + mfj + 
m 2 H , hence the effect of RGE running is the same as in mSUGRA with the replacement 
m 2 — > m 2 . If sin (3 ~ 1, correct gauge symmetry breaking requires /i 2 ~ —vr? Hu — M%/2, where 
all quantities are taken at the weak scale. Hence 5m 2 H directly affects the value of at the 
weak scale. For the string scenario (§) the effect is limited because the total scalar mass is 
bounded from above by the gaugino mass, but for the 50(10) scenario @ the change of |/z| 
from its mSUGRA prediction can be more drastic. 

Fig. 3 shows the correlation between Aa^ and a xv for tan/3 = 10 and fi > in three 
different models: mSUGRA, the string scenario described by eq.(D) with universal sfermion 
masses (= m), and the 50(10) GUT model. Note that the assumption of universal sfermion 
masses implies rn? H = rn 2 Hu at the GUT scale in the stringy model; the value of the Higgs 
soft breaking masses is determined by eq.([3j) for given m and M, and A = —M. For the 
5O(10) GUT model, we assume m 16 = m 10 = m and fix A = and M|, = —50000 GeV 2 . 
As before, we scan the region m < 500 GeV and M < 600 GeV, requiring the physical Higgs 
mass and chargino mass to lie above their experimental lower bounds, and excluding regions 
with charged LSP. The heavily marked points again indicate the cosmologically favored region 
where Vt x h 2 < 0.3. 

In our superstring scenario the higgsino component of the LSP can be larger or smaller than 
in mSUGRA, depending on the ratio m/M. For fixed M, a small m gives a larger higgsino 
component of the LSP because then m# should be larger than m 2 in order to satisfy eq.@, 
i.e. 5m 2 H > in eq.(||]). This in turn reduces the supersymmetric contribution fi 2 to the Higgs 
mass term which is required to obtain the correct Z boson mass. On the other hand, for 
m 2 > M 2 /3, eq.(H|) gives 5m 2 H < 0, leading to even larger than in mSUGRA. Altogether 
this model allows for a broader range of values of a xp , with maximum around 1.7 x 10~ 8 pb, 
about a factor of 2 above the maximal mSUGRA prediction. On the other hand, the maximum 
value of Aa M is reduced to 2.3 x 10~ 9 from 3.3 x 10 -9 in mSUGRA. This reduction mainly 
comes from the fact that the minimal allowed value of M is raised by about 40 GeV due to 
the choice A = —M (rather than A = in the mSUGRA case); we saw in Fig. 2 that negative 
values of A tend to reduce if all other input parameters are kept fixed. The variation of \x 
is less important for Act^. The reason is that increasing \\x\ increases the higgsino mass and 

7 The effective en is slightly negative for low SUSY breaking scale, but turns positive if this scale is large. 
Recall that the relevant scale for the analysis of gauge symmetry breaking increases with increasing sparticlc 
masses. 
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Figure 3: Aa M vs. a xv in three different models : mSUGRA, a superstring scenario and the minimal SOiVS) 
GUT model. We chose tan/? = 10 and /i > 0. The scanning procedure and notation is as in Fig. 1. 



reduces higgsino-gaugino mixing, which reduces |Aa M |, but at the same time increases [Il — ^r 
mixing, which tends to increase |Aa M |. 

We also performed a scan (not shown) for a second ansatz for scalar soft breaking masses 
which is compatible with the general expressions (||). Here we kept the soft breaking con- 
tributions to the squared squark and Higgs boson masses equal to each other. This implies 
m Q = m u = m D = m H u = m H d = M 2 /3 and m 2 E + m\ = 2M 2 /3. We then allowed vtie 
and mi to differ from each other, keeping universality among different generations. This scan 
gives results for a xp which are quite similar to those of the mSUGRA scenario; however, the 
predictions for a M are generally lower than in mSUGRA. This is true in particular for param- 
eter values that satisfy Q x h 2 < 0.3. For the given relation between average slepton mass and 
gaugino mass, an acceptable value of the relic density can only be achieved^] if m E < 
note that due to its larger hypercharge, exchange contributes 8 times more strongly to the 
annihilation of bino-like LSPs than plus v exchange does. This mass splitting between <3r 
and v reduces Aa M , since the (usually dominant) sneutrino-chargino loop contributions are 
suppressed compared to the scenario with universal slepton masses. The main difference to 
our first stringy scenario, where we allowed different soft breaking contributions for sfermions 
and Higgs bosons while keeping sfermion masses equal to each other, are the lower maximally 
allowed values of both a xp and Aa M . This also leads to smaller values of Aa M and a xp in the 
cosmologically acceptable region: the minimal a xp which is compatible with Q x h 2 < 0.3 is 
reduced to ~ 2.5 x 10~ 9 pb from ~ 6 x 10~ 9 pb, while the smallest cosmologically allowed 
Aa M is reduced from ~ 1.25 x 10~ 9 to ~ 0.9 x 10~ 9 . 

While the predictions of these versions of the superstring inspired model differ from those 
of mSUGRA "only" by about a factor of 2, the 5O(10) GUT model with M 2 D = -50000 GeV 2 

8 We note in passing that the dilaton-dominated scenario pq , where all scalar soft breaking masses are 
equal to M/y/3, cannot simultaneously accommodate > 111 GeV and Q x h 2 < 0.3, unless tan/3 is very 
large. 
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leads to a dramatically changed correlation between Aa M and a xp . Since rn 2 Hu at the GUT 
scale is larger than in mSUGRA, the higgsino component of the LSP increases. This effect is 
especially significant for small M, where the gaugino mass contributions to the RGE of the 
scalar masses are not so large, so that a modification of the GUT scale boundary condition 
more strongly affects weak scale values of scalar masses; see eq.(|). In this region a XP increases 
quite a lot, reaching values up to ~ 4 x 10~ 7 pb, which could be detected in the near future. 
In contrast, the allowed values of Aa M are less than 1 x 10~ 9 , well below near future sensitivity. 
This reduction is caused by the increase of the sneutrino mass and the decrease of the jxn mass 
due to the D— term contributions in eq.([|). Because of this large slepton mass splitting, eq.(|TJ) 
is not applicable in this scenario. The D— terms suppress the contribution from chargino- 
sneutrino loops, while enhancing the contribution from neutralino-smuon loops. These two 
contributions can now have comparable magnitude, but opposite sign, leading to significant 
cancellations in the small m region. These cancellations explain why we sometimes find an 
anti-correlation between a xp and Aa^ in this model, i.e. decreasing m 2 can lead to decreasing 
Aa M (and increasing a xp , as usual). Of course, for large m all contributions to Aa M are small. 



<3 



°xp (pb) 



Figure 4: Aa M vs. a xv in mSUGRA and the 50(10) model with Mf, = -50000 GeV 2 . Here, we fix m and 
M such that m x o ~ 120 GeV and m^ R ~ 182 GeV, and vary A and tan j3. 

In Fig. 3 the models with and without SO (10) D— terms are well separated. Therefore 
one might conclude that one can test for the existence of 5*0(10) D— terms by measuring Aa M 
and/or a xp . However, in Fig. 3 we fixed the values of tan/3, A and Mp. In order to see the 
effects of varying tan/3 and A, we fix m = 140 GeV and M = 280 GeV in mSUGRA, and 
m = 260 GeV and M = 300 GeV in the 5O(10) model, so that m^o ~ 120 GeV and m- eR ~ 182 
GeV in both cases. We then scan the region of — 2M < A < 2M and the entire allowed range 
of tan/3. Both the lower and upper bounds on tan/3 come from Higgs searches at LEP; recall 
that in models with radiative symmetry breaking, all neutral Higgs bosons become quite light 



as tan/3 approaches its upper bound |28|. Fig. 4. shows the result of this scan. Here the 



short-dashed line corresponds to A = 0, with positive (negative) values of A falling to the left 



10 



(right) side of this line. As expected, Aa M increases when tan/3 increases, and depends only 
weakly on A. On the other hand, a xp increases with increasing tan/3, but decrease when A is 
increased. As noted earlier, this reduction of o XP comes from the fact that the Higgs masses 
increase, and the LSP becomes more bino-like, when A increases. 

In this figure, mSUGRA scenarios are characterized by relatively large Aa M and small a xp , 
but the opposite is true for the SO (10) GUT model with M 2 D = -50000 GeV 2 in most of the 
allowed region. The two classes of models thus remain quite well separated. Of course, it 
is by no means certain that the extra D— term contribution in the 5*0(10) model should be 
large and negative. In fact, if one requires that all three third generation Yukawa couplings 
unify at the GUT scale, which is not unreasonable in minimal 5*0(10) since all third generation 
superfields reside in a single representation of 5O(10), consistent radiative symmetry breaking 
is only possible if M|> > H38| . In Figs. 5a, b we therefore show the dependence of Aa M and 
a xv on Ml, for m = 260 GeV, M = 300 GeV, A = and tan/3 = 10. The upper and 
lower bounds on come from the requirement that SU(2) singlet and doublet sleptons, 
respectively, should be heavier than the lightest neutralino. We see that a xp increases very 
quickly for negative M|,; near the lower bound on Mp one even approaches the sensitivity 
of current Dark Matter search experiments. Conversely, Aa M increases if M|> is raised above 
0, but the maximal increase amounts to less than a factor of 2 compared to the mSUGRA 
situation, which corresponds to M|> = 0. Note, however, that Yukawa unification also requires 
tan/3 > 35 [38|; the combined enhancement of Aa^ from large tan/3 and from M|, > implies 
that measurements of the anomalous dipole moment of the muon might be one of the best 
ways to test 5O(10) models with Yukawa unification. 




a) b) 

Figure 5: The dependence of a) a xp and b) Aa^ on M% in the 5*0(10) model, for m = 260 GeV, M = 300 
GeV, A = and tan/3 = 10. 
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4 Conclusions 



In this paper, we investigated the SUSY contribution to the muon MDM a M , and to the 
spin-independent neutralino-proton cross section a xp , in several SUSY models with either 
universal or non-universal soft scalar masses at the GUT scale. Both these observables can 
become significant if tan/3 is large, but will be difficult to measure in the small tan/3 region. 
Both quantities are sensitive to chirality violation in the matter (s)fermion sector, which in 
the MSSM is enhanced for large tan/5. Conversely, if tan/3 is small the experimental Higgs 
mass bounds tightly restrict the allowed SUSY parameter space, forcing the SUSY breaking 
scale to be quite high. While these general statements are fairly model-independent as long as 
one sticks to the field content of the MSSM, quantitative predictions do depend significantly 
on details of the spectrum of superparticles, in particular on the implementation of radiative 
gauge symmetry breaking. 

We considered three different models: mSUGRA, a superstring scenario where SUSY is 
broken by the F— terms of the dilaton and the moduli superfields, and an 5*0(10) GUT model. 
The superstring model in principle has a very large parameter space, since it allows for non- 
universal masses for matter sfermions and Higgs bosons, subject only to the sum rule (|2|). Here 
we explored two orthogonal directions in this parameter space, where universality is violated 
either only in the Higgs sector or only in the slepton sector. We found that allowing the Higgs 
masses to differ from the sfermion masses at the GUT scale can vary the predicted value of 
o XP by a factor of a few in either direction, compared to the mSUGRA prediction for the 
same parameters. Given the uncertainty in the prediction of a xp due to unknown hadronic 
matrix elements, and the uncertainty in the predicted LSP-nucleus scattering rate due to the 
uncertainty of the ambient LSP flux, such a variation is barely significant. The prediction of 
Aa M is almost insensitive to this variation of the Higgs soft breaking masses. Conversely, if we 
maintain the universality of the soft breaking contributions to squark and Higgs masses but 
allow the masses of SU(2) doublet and singlet sleptons to differ, we can change the prediction 
for Aa M by a factor of ~ 2 around the mSUGRA value, while keeping a xp constant; the 
predicted LSP relic density, and hence the cosmologically preferred region of parameter space, 
also depends on this ratio of slepton masses. Generally speaking, in mSUGRA as well as its 
stringy variant the measurement of the anomalous dipole moment of the muon seems to hold 
more promise in the near future than searches for LSP Dark Matter; this is largely due to the 
most recent severe constraints on the Higgs sector from LEP experiments. 

We should caution the reader that this conclusion crucially depends on our choice to 
only study scenarios where both squarks and gluinos can be produced at the LHC. Given 
current slepton and chargino search limits, this implies that sfermion and gaugino masses are 
very roughly of the same order of magnitude. If one allows scalar masses m > 1 TeV and 
m 2 ^> M 2 , one can find solutions in mSUGRA where the LSP has a large, or even dominant, 



higgsino component [39]. Since in this region of parameter space all sfermions are very heavy, 



the supersymmetric contribution to a M is totally negligible, but a xp can be very large ||40| 



Much more dramatic deviations from mSUGRA predictions are possible if one introduces 
SO(10) D— term contributions to scalar masses. In this model the maximal allowed value of 
a xp for given tan/3 can exceed the mSUGRA prediction by more than a factor of 30, if the 
D— term contribution to the mass of the Higgs boson that couples to the top quark is positive. 
For the same sign of the D— terms the supersymmetric contribution to is reduced by a 
factor of a few, since the change of the slepton masses increases the importance of neutralino 
loop diagrams compared to chargino loops, leading to a strong cancellation in the total result. 
In this kind of model direct Dark Matter searches therefore appear to be more promising in 
the next few years than measurements of a M . 



12 



A striking result of our analysis is that even if we restrict squark and gluino masses to lie 
within discovery range of the LHC, and in addition fix the value of tan/3, the prediction for 
o XP still varies by almost three orders of magnitude within the range of models we studied; in 
contrast, the prediction for Aa M "only" varies by about one order of magnitude; see Fig. 3. 
One reason for this difference is that a xp measures squared Feynman amplitudes, while is 
directly proportional to a (loop) amplitude. In addition, we saw that Aa^ does not depend on 
\i very sensitively, since various effects tend to cancel. In contrast, the sensitivity of o XP to fi 
is often enhanced by cancellations between diagrams with different fi— dependence. A similar 
remark holds for the tan/? dependence. Fig. 1 shows that it is "only" linear for Aa^, but much 
stronger for a xp . In this case the strong sensitivity of a xp is partly due to the dependence of the 
mass of the heavy CP-even Higgs boson on tan/3, which is a direct consequence of radiative 
gauge symmetry breaking. Given the theoretical uncertainties involved in predicting Dark 
Matter detection rates, this strong parameter dependence of a xp is in some sense fortunate; 
it means that even an order-of-magnitude determination of a xp will allow us to make signifi- 
cant statements about supersymmetric parameters. Given its comparatively mild parameter 
dependence, it is also fortunate that the prediction of is much cleaner, if we assume that 
the hadronic uncertainty can be reduced to the level of the expected experimental uncertainty, 
which seems feasible. 

Our overall conclusion is that measurements of and/or a xp might well yield a positive 
signal in the next few years. These measurements can help to distinguish between different 
supersymmetric models. Since well motivated models exist where one of these quantities is 
sizable while the other is small, it is important to measure, or at least constrain, both of 
them as precisely as possible. On the other hand, we also saw that two models with very 
different scalar spectrum (the 50(10) model and mSUGRA with m 2 3> M 2 ) can lead to very 
similar predictions for both a M and o XP . This once again illustrates the truism that indirect 
measurements can provide crucial information, but they can never completely replace direct 
searches for new particles. 
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